Weak ferromagnetism with very large canting in a chiral lattice: (pyrimidine) 2 FeCl2 
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The transition metal coordination compound (pyrimidine) 2 FeCh crystallizes in a chiral lattice, 
space group 74i22 (or 74322). Combined magnetization, Mossbauer spectroscopy and powder neu- 
tron diffraction studies reveal that it is a canted antiferromagnet below Tjv = 6.4 K with an unusually 
large canting of the magnetic moments of 14° from their general antiferromagnetic alignment, one of 
the largest reported to date. This results in weak ferromagnetism with a ferromagnetic component 
of 1 [is- The large canting is due to the interplay between the antiferromagnetic exchange inter- 
action and the local single-ion anisotropy in the chiral lattice. The magnetically ordered structure 
of (pyrimidine)2FeCl2, however, is not chiral. The implications of these findings for the search of 
molecule based materials exhibiting chiral magnetic ordering is discussed. 

PACS numbers: 75.25.+z; 75.50.Ee; 75.50.Xx 



I. INTRODUCTION 

One major research effort in the very active field of 
molecule based magnetic material9i*2i2i2£ is directed to- 
wards synthesizing multifunctional compounds that com- 
bine magnetism with a second physical property, such as 
conductivity^ or optical activity^ In this framework, ex- 
tensive efforts are undertaken to find chiral magnetic ma- 
terials that exhibit magnetochiral dichroism (MChD), a 
phenomenon first observed by Rikken and Raupach in a 
chiral paramagnetic material* 8 Large MChD is expected 
in materials that combine chirality and magnetic order. 
To date, however, there is no evidence for large MChD, 
although a number of chiral magnetic materials have been 
report ediS 

Recently, the magnetism of pyrimidine bridged tran- 
sition metal complexes has been investigated in which 
pyrimidine (C4H4N2), in the following abbreviated 
as PM, plays the role of an antiferromagnetic cou- 
pling unitii2iiiii2ii£ Interestingly, the halide complexes 
(PM) 2 TX 2 (T = Co 11 , X = CI, Br; T = Fe 11 , X = CI) 
possess a chiral 3-dimensional network of T ions and ex- 
hibit weak ferromagnetism below about 5 Kii^i^ The 
preliminary analysis of magnetization measurements on 
these compounds pointed to a very large canting and an 
associated large ferromagnetic component of the ordered 
moments. It appeared possible that these compounds 
possess a chiral magnetically ordered structure. 

In order to determine the magnetic structure and to 
elucidate the origin of the large canting we performed 
magnetization measurements, Mossbauer spectroscopy 



and powder neutron diffraction on (PM)2FeCl2. Here, we 
present a complete analysis of these data together with 
electronic structure calculations. We discuss the impli- 
cations of our findings for studies in search of molecule 
based materials exhibiting chiral magnetic ordering. 



II. EXPERIMENTAL 

Microcrystalline samples of (PM)2FeCl2, in the follow- 
ing abbreviated PFC, were obtained by mixing aqueous 
solutions of Fe 11 chloride (e.g. 10 mmol in 20 ml H 2 0) 
with the stoichiometric amount of pyrimidine. The re- 
sulting yellow precipitates were filtered and washed thor- 
oughly with H 2 0. Mossbauer data indicate the presence 
of a secondary phase, which we identified as (PM)FeCi2. 
The volume amount of the secondary phase increases 
with storage time of the samples. While a newly made 
sample contained less than 2% secondary phase, about 
22% secondary phase have been observed in the same 
sample stored for two years under air. This indicates 
that under such conditions PFC is unstable against the 
formation of (PM)FeCi2. Magnetization measurements 
were carried out using an MPMS Squid magnetometer 
(Quantum Design). Some 20 mg of freshly prepared sam- 
ple was filled into a gelatine capsule. 57 Fe Mossbauer 
spectroscopy experiments have been executed in a stan- 
dard low-temperature Mossbauer set-up at temperatures 
ranging from 3 K to 300 K (source: 57 Co-in-Rh matrix at 
room temperature; half width half maximum (HWHM): 
0.130(2) mm/s). The spectra have been evaluated us- 



ing the Mossbauer fitting program Recoil^ in the thin 
absorber approximation. Above Tjv the spectra were 
modeled as Lorentzian lines in the presence of an elec- 
tric field gradient. Below Tjv we employed a method for 
the calculation of Mossbauer line shapes in the presence 
of magnetic dipole and electric quadrupole hyperfine in- 
teractions introduced by Blaes et alH Neutron powder 
diffraction measurements were performed using the in- 
struments E6 and E9 at the Berlin Neutron Scattering 
Center (BENSC). Neutron wavelengths of 2.448 A and 
1.7964 A, respectively, were used. The instrument E6 
provides a high neutron flux and medium resolution, is 
equipped with a 20°-multichannel detector, and covers a 
range of scattering angles up to about 100°. In contrast, 
the instrument E9 is a low-flux high-resolution powder 
diffractometer with an extended 26>-range up to 160°. 
Therefore, the former was used for the study of mag- 
netic Bragg reflections, while the latter was employed for 
checking the crystal structure at low temperatures. The 
non-deuterated sample was filled into a 7 mm diameter 
vanadium can with 40 mm length, resulting in a sample 
volume of about 1.5 cm 3 . An absorption correction for 
cylindrical samples was applied (fiR = 0.94) to account 
for the strong incoherent scattering from hydrogen. Dur- 
ing the refinement, no other but symmetry constraints 
were used for the atom positional parameters. There- 
fore, the validity of the crystal structure model is proven 
by the correct geometry of the pyrimidine molecules de- 
termined in the refinement. The Rietveld refinement 
of the diffraction data was carried out using the WIN- 
PLOTR/FULLPROF packaged 
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FIG. 1: Crystal and magnetic structure of PFC as determined 
from low-temperature powder neutron diffraction. Lower sec- 
tion: view along the c axis with the positions of the 90° screw 
axes marked. The angle 59° between the magnetic moments 
and the local EFG is also marked. 



IV. RESULTS AND DISCUSSION 



III. THEORETICAL CALCULATIONS 

The electronic structure calculations were performed 
with a code based on a local basis set^S The unrestricted 
Hartree-Fock method and the hybrid functional B3LYP 
is applied for solids, and the basis functions are chosen as 
Gaussian type orbitals. The iron basis sei£& (outermost 
d-exponent 0.4345) is of the size [5s4p2d], the chlorine 
basis seliSi (with one (i-exponent with value 0.5) of the 
size [5s4pld], carbon and nitrogen basis sets^ of the size 
[3s2pld], and finally a [2slp] hydrogen basis selpS 3 - was 
used. To test the stability of the results, additional tight 
basis functions were added for the iron atoms and a dif- 
fuse sp-function (exponent 0.12) was added at the nitro- 
gen site to account for a better description of this neg- 
atively charged atom. The results were, however, found 
to be essentially stable with respect to the various basis 
sets. With these parameters, the Hartree-Fock or Kolm- 
Sham equations for a ferro- or antiferromagnetic (Neel- 
like) structure are solved self-consistently, and properties 
such as charge distributions and field gradients can be 
computed. The structural data used in the calculations 
is based on the measurements as described in the follow- 
ing. 



A. Crystal structure 

The high-resolution powder neutron diffraction data, 
taken at 10 K, confirm that PFC is isostructural to 
the Co analogue^, namely tetragonal with the chiral 
space group 74i22 (or the enantiomorphic 74322) with 
a = b = 7.4292(3) A, c = 20.364(1) A. The Fe-N and 
Fe-Cl distances 2.26 and 2.39 A are determined, respec- 
tively. The crystal structure is shown in Fig. 1. It con- 
sists of a chiral three-dimensional network of Fe ions co- 
ordinated by two CI and linked by pyrimidine. Due to 
the 4i (or 43) screw axis, the local environment of two 
Fe ions in neighboring layers z and z + 1/4 is rotated 
by 90°. We will discuss below that this specific feature 
of the crystal structure of PFC can be regarded as the 
origin of the large canting observed in the magnetically 
ordered state. 

In the local FeN4Cl2 geometry all nitrogen atoms of 
pyrimidine are equatorially coordinated. It was pointed 
oulii^ that the electron spins in magnetic d x 2_ y 2 or- 
bitals are antiferromagnetically correlated through the 
pyrimidine molecular orbital(s), i.e., by a superexchangc 
mechanism (SiS^LSi when the nitrogen lone pair and 
d x 2_ y 2 orbitals have an appreciable orbital overlap on 
both sides. 
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V. RESULTS AND DISCUSSION 



A. Magnetization 

Figure 2 shows the M vs. H hysteresis curve of PFC 
measured at 2.0 K. For clarity, the virgin curve is not 
shown. On increasing the field from zero, we observe an 
initial steep increase of the magnetization in a low field 
region (a few 100 Oe) to a value much smaller than the 
full Fe 11 moment, and subsequently a rounded crossover 
around 500 Oe to a linear field dependence up to 55 kOe 
(not shown). Such a behavior is typical of a weak fer- 
romagnet, where the steep increase at low T reflects the 
spontaneous magnetization and the linear high-field be- 
havior is due to strong antiferromagnetic interactions. 
On sweeping the field down from 55 kOe, a weak hys- 
teresis develops and a sharp kink is observed at H = 0. 
This kink allows an accurate determination of the spon- 
taneous magnetization. We determine a value of 0.31(1) 
/jg on the powder average. The coercitive field is small, 
150 Oe, and therefore PFC behaves as a soft magnet. 
Assuming that the ferromagnetic component in a given 
magnetic domain is confined to a specific crystal direc- 
tion, the measured powder average of the spontaneous 
magnetization has to be multiplied by a factor of three to 
obtain the ferromagnetic moment along this axis. This 
fact is frequently overlooked by other authors. Thus, 
the magnetization measurements result in a ferromag- 
netic component of 0.93(3) /is per Fe ion in the canted 
antiferromagnetic state. We will show in the following 
that this value is in full agreement with the combined 
Mossbauer spectroscopy and neutron diffraction results. 



B. Mossbauer spectroscopy 

In Fig. 3 we plot the Mossbauer-spectra of PFC at 
temperatures between 3.75 K and 300 K. At highest tem- 
peratures we observe a two line spectrum, resulting from 
a quadrupole splitting of QS — 3.112(3) mm/s and an 
isomer shift of IS = 0.909(1) mm/s (relative to the 57 Co- 
in-Rh source). These values are typical for high spin Fe 11 . 
Upon lowering the temperature, and above Tjy, both QS 
and IS slightly increase, and we find at 10.4 K values of 
QS = 3.271(2) mm/s and IS = 1.037(1) mm/s. 

Below T^, at 3.75 K, eight separate absorption lines 
are observed. These we attribute to the coexistence of an 
electrical field gradient (EFG) and a hyperfine magnetic 
field (£?hf)- A fit to the data, taking into account the 
EFG, -Bhf and the asymmetry parameter rj = (V xx — 
V yy )/V zz yields values of IS = 1.035(3) mm/s, r\ — 0.106, 
QS = 3.28(1) mm/s and a hyperfine magnetic field of 
30.19(3) T as best solution. Moreover, from this fit we 
determine a temperature independent angle of -Bhf with 
the main EFG component V zz of 59.3(1)° at each Fe site. 




H(kOe) 

FIG. 2: Low-field section of the hysteresis-loop of PFC mea- 
sured at 2 K. The field was swept between -55 kOe and 55 
kOe. For clarity, the virgin curve is not shown. Horizontal 
arrows mark sharp kinks allowing for the determination of the 
spontaneous magnetization. 

C. Results from calculations 

From the calculations, the charges were determined to 
be +1.9 for Fe, -1.0 for CI, -0.8 for N; between and 
+0.7 for carbon, and the hydrogen atoms were found to 
be slightly positive charged (at most +0.1). The spin at 
the Fe sites is 1.9 and the magnetic moment 4.04 fig, us- 
ing a ^-factor of 2.13^ Thus, there is moment reduction 
by about 5% due to covalency effects. The magnitude 
of the spin at the CI, N or C sites is less than 0.1 and 
thus negligible, and the spin at the hydrogen sites virtu- 
ally zero. The electric field gradient at the Fe site has 
components with the value V xx — —0.8, V vv = —1.3, and 
V zz = 2.1 x 10 22 V/m 2 . These numbers have an uncer- 
tainty of the order of 0.1 x 10 22 V/m 2 . The asymmetry 
parameter rj is especially sensitive to this uncertainty and 
values in the range from 0.1 to 0.3 are obtained. 

The second largest component of the tensor is identical 
with the crystallographic c-direction, the other two com- 
ponents lie perpendicular in the tetragonal basal plane 
and are rotated by 45° with respect to a and a'. The 
largest component is in the direction of the CI atoms. 
We can compare the computed field gradient with the 
experimental value by converting according to the for- 
mula QS = eQV zz c/2E , with E = 14.4 keV and using 
the value of 0.16 barn for the quadrupole moment Q. 28 
With these parameters, we obtain a computed QS = 3.5 
mm/s which is in good agreement with the experimen- 
tal value. Finally, it is worthwhile mentioning that the 
aforementioned calculated values hardly depend on the 
magnetic ground state, and similar values are obtained if 
a ferromagnetic or an antiferromagnetic ground state is 
assumed. 

An attempt to extract information about the magnetic 
exchange coupling could be made by comparing energies 
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FIG. 3: Mossbauer spectra of PFC measured at various tem- 
peratures. Arrows mark eight absorption lines observed at 
3.75 K . The solid lines are fits to the data points (see text). 



for ferro- and antiferromagnetic structure. This differ- 
ence is (per Fe ion) roughly 0.01 eV at the HF level (with 
the ferromagnet being lower), and roughly 0.08 eV at 
the B3LYP level (with the antiferromagnet being lower). 
However, these energies should be of the order of the 
Neel temperature (6.4 K or 0.6 meV). The large devia- 
tion could be due to the fact that a simplified magnetic 
order is used in the simulations, neglecting the canting, 
or that the description of the electronic correlations is too 
crude due to the approximations in the functional. If the 
HF result was assumed to be reasonable, then the effect 
of electronic correlations would have to be of virtually the 
same size, but with opposite sign. Such a subtle effect 
is very difficult to be computed properly, and a density 
functional calculation may be too crude. 



D. Magnetic structure 

Figure 4 shows the difference of the powder neutron 
diffractograms taken at 1.6 and 7 K (i.e., well below 
and above Tjy), revealing the additional Bragg intensi- 
ties produced by the long-range magnetic ordering. All 
magnetic reflections can be indexed on the basis of the 
crystallographic unit cell. The rule h + k + I = even for 
the observed magnetic Bragg reflections (hkl) indicates 
that the body centering is conserved. The presence and 



strengths of the reflections (00Z) suggest that the mag- 
netic moments lie basically perpendicular to the c axis 
and that moments in neighboring layers z and £+1/4 are 
basically oriented antiparallel. Therefore, the antiferro- 
magnetic coupling through pyrimidine bridges is clearly 
indicated. 

First refinements of the neutron diffraction data re- 
sulted in an ordered moment of roughly 4 /is at 1.6 K. 
However, the orientation of the moments within the basal 
plane as well as the angle and direction of canting can not 
be determined from these data since the only clear sig- 
nature of the ferromagnetic component is the weak (112) 
reflection which is hidden in the noise. We therefore used 
the information from the Mossbauer spectroscopy to ar- 
rive at the final magnetic structure model. 

Assuming the electric field gradient parallel to the Fe- 
Cl bonds, as it was determined in the calculations, in 
the basic antiferromagnetic structure the moments have 
to be aligned parallel to the a axis in order to enclose 
an identical angle with the local EFG at each Fe. To 
increase this angle from 45° to 59.3° - as observed in the 
Mossbauer experiment - either a canting by 14.3° within 
the basal plane (along a') or by 45° along c has to be 
introduced. However, only a canting in the basal plane 
gives a ferromagnetic component close to 1 [Lb , whereas a 
canting along c would require a ferromagnetic component 
of 2.8 hb- Comparing these values to the results from 
the magnetization measurements, the latter model can 
be clearly ruled out. Therefore, in the final magnetic 
structure model, the moments were confined to the basal 
plane and the canting angle 14.3° was fixed. The solid 
line in Figure 4 is a fit of this final model to the data. 
The only free parameter is the magnitude of the ordered 
moment. All other parameters, such as lattice constants 
and lineshapes, were taken from the fit of the nuclear 
Bragg pattern recorded at 7 K (not shown) . 

Given the large background noise from the hydrogen 
incoherent scattering, it gives a good agreement with the 
data (R p = 0.114, R wp = 0.137). The fit yields an or- 
dered moment of 4.0(3) [Ib - in agreement with the above 
calculated full moment on the Fe - and a ferromagnetic 
component of 1.0(1) [is- The latter value is in good 
agreement with the magnetization data. 

The temperature dependence of the intensity of the 
(002) Bragg reflection is shown in Figure 5 together 
with the corresponding data, the square of the hyper- 
fine field -Bhf, from the Mossbauer experiment. It de- 
creases continuously with rising temperature and van- 
ishes at 6.4(2) K. Assuming that the canting angle is con- 
stant, the Bragg intensity is proportional to the square 
of the ordered moment. A fit of the region close to Tjv 
with a power law is consistent with a critical exponent of 
(3 = 0.36(1) and therefore with 3D Heisenberg behavior. 

The canting angle a = 59.3° - 45° = 14.3° observed 
in the present measurements is extremely large, actually 
one of the largest reported to date for any weak ferro- 
magnetic compound. For example, weak ferromagnetism 
of Fe with a canting angle of 16° has been reported for 
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the intermetallic compound UFe4Al§. However, the fer- 
romagnetic component is only 0.3 \xb per Fe in this com- 
pound and the canting is due to an interaction between 
the U 5/ and the Fe 3d electrons^ Therefore, the physics 
of this compound is hardly comparable to that of PFC. 

Large canting angles of 2-7° have been also observed in 
Gd 2 Cu04-type cuprates. For these materials, a correla- 
tion between weak ferromagnetism and the crystal sym- 
metry has been discussed in detail^ We believe that 
such a correlation is also present in PFC and that the 
large canting observed in this compound is a direct re- 
sult of its chiral lattice symmetry. Due to this sym- 
metry, the orientation of the local easy axis varies by 
90° between nearest Fe neighbors linked by a pyrimidine 
molecule. Therefore, the local anisotropy would favor 
a 90° alignment between the moments. However, the 
AFM interaction favors a 180° alignment. The actual 
angle 180° — 2a = 151.4° between neighboring moments 
therefore can be regarded as the result of the compe- 
tition of the local anisotropy and the antiferromagnetic 
exchange. The final magnetic structure model suggests 
that the local easy axis lies perpendicular to c and within 
the equatorial plane of the local coordination octahedron 
of each Fe. Thus, it coincides with the axis of the smallest 
EFG component V xx calculated above. We may express 
this anisotropy as an additional term in the Hamiltonian 
of the system, 



U = J2[JSi-S 2 +D{S lx S 2y ) 



(1) 



NN 



where the sum is over all pairs of nearest neighbor (NN) 
spins (SijS^), Sj = (Sj X ,Sjy,Sj Z ), and D < 0. The re- 
sulting angle a is related to D and J by D/ J — —2 sin 2a. 
With the measured a we get D = — 0.96J. If one would 
assign the canting to a Dzyaloshinsky-Moriya (DM) type 
interactionj2ii22i2ii the DM vector was aligned along c as 
expected for symmetry reasons and the DM term would 
read D(S\ X S2 V — S\ y S2 X )- The latter term of (1) differs 
from the DM term in leading to a canting towards spe- 
cific axes x and y for Si and 5*2, respectively, rather than 
a canting within the tetragonal basal plane in general. 

It is interesting to note that magnetization data on the 
Co analogues to PFC, (PM) 2 CoCl 2 and (PM) 2 CoBr 2 , 14 
point to even larger ferromagnetic components and there- 
fore even larger canting angles than in PFC. Another ex- 
ample of such a large canting may be 3D-[Fe(N3) 2 (4,4'- 
bpy)], with bpy = 4,4'-bipyridine. A microcrystalline 
sample of this compound was recently reported to ex- 
hibit a large spontaneous magnetization of 0.48 hb per 
Fe tentatively ascribed to ferromagnetic ordering while 
a canted structure was not excludedi^i^ This tetrag- 
onal compound also forms a chiral 3D network struc- 
ture, space group P4i2i2, which is closely related to 
that of (PM) 2 FeCl 2 . The negative Curie- Weiss temper- 
ature and the similarity of the magnetization data of 
3D-[Fe(N 3 ) 2 (4,4'-bpy)] with that of the (PM) 2 TX 2 com- 
plexes suggests that also in the former a canted antifer- 
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FIG. 4: Difference of the two neutron powder diffractograms 
measured at 1.6 K (below T N ) and 7.0 K (above T N ) for PFC. 
The wavelength was 2.448 A. The reflections are indexed on 
the basis of the crystallographic unit cell. 
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FIG. 5: Temperature dependence of integrated intensity of 
the (002) Bragg reflection, being proportional to the square 
of the ordered moment, together with the hyperfine field Bhf 
observed by Mossbauer spectroscopy. The dashed line is a 
guide to the eye. The solid line is a fit to a power law which 
is consistent with 3D Heisenberg behavior. 



romagnetic state with a very large canting (of roughly 
20°) within the tetragonal basal plane is realized. 

To our knowledge, only few other molecule based 
magnets with chiral crystal structures have been re- 
ported. The oxalato (ox) based compound [Co(2,2'- 
bpy)s] [Co 2 (ox)3]C104^i exhibits weak ferromagnetism 
with small canting (/if ro = 0.009 (J.b)- The compounds 
of the series [Z n (2,2'-bpy) 3 ] [C10 4 ] [M n Cr m (ox) 3 ] 37 order 
ferromagnetically. Both systems crystallize in the cubic 
chiral space group P4i32. Another example is ferrimag- 
netic Mn(hfac) 2 NITPhOMe, crystallizing in space group 
P3 1 $L The other compounds discussed in the frame- 
work of chiral magnetism are based on chiral constituents 
and crystallize in non-centrosymmetric but achiral space 
groups, such as PI or P2\2\2\? 
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VI. CONCLUSION 

We have shown that (PM)2FeCl2 is a weak ferromagnet 
with low coercitive field but a very large ferromagnetic 
component of 1 fig, corresponding to a canting angle of 
14°. This is one of the largest values reported to date. 
In a given magnetic domain, the antiferromagnetic com- 
ponent of the moments is confined to the a axis, while 
the canting occurs along the perpendicular direction a' . 
We argued that the large canting in PFC is the direct 
result of its chiral lattice symmetry leading to a compe- 
tition between the antiferromagnetic exchange interac- 
tion and the local single-ion anisotropy. The appearance 
of canting rather than a chiral magnetic structure may 
be a general feature of the magnetic ordering in simi- 
lar structures possessing 90° screw axes. This observa- 
tion is of interest in view of the extensive efforts to pro- 
duce a molecular magnet exhibiting chiral magnetic or- 
dering. To date, both strategies, crystallizing (i) achiral 
constituents in a chiral lattice and (ii) chiral constituents 



in non-centrosymmetric but achiral lattices did not lead 
to chiral magnetic ordering. 

Most of these compounds possess 90° or 180° screw 
axes. These lead to a ±90° or to no alteration at all, 
respectively, of the local anisotropy along the screw axes 
and thus do not necessarily support any chirality of the 
ordered magnetic structure. Local anisotropy supporting 
chiral magnetic ordering therefore appears more likely 
in chiral trigonal/hexagonal lattices. We suggest that 
the search for chiral magnetic ordering should focus on 
compounds of that kind. 
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